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Abstract. We mimic the stochastic Hamiltonian reduction of Lazaro-Cami and Ortega [17, 18] for the case 
of certain non-holonomic systems with symmetries. 

Using the non-holonomic connection it is shown that the drift of the stochastically perturbed n-dimcnsional 
Chaplygin ball is a certain gradient of the density of the preserved measure of the deterministic system. 
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1. Introduction 

Imagine a ball sitting on a rough horizontal table. Because of the roughness of the table this ball -the 
Chaplygin ball- cannot slip, but it can turn about the vertical axis without violating the constraints. Its 
geometric and gravitational center coincide but there may be an inhomogeneous mass distribution. Suppose 
now that the ball is subjected to a Brownian noise such that there is a random jiggling in all of its angular 
and translational degrees of freedom. This problem is similar to the stochastic rigid body considered in [18] 
but upon imposing the no-slip constraints some differences can be expected. One may wonder if the stochastic 
Chaplygin ball will acquire a drift that makes it roll on the table or spin about its vertical axis or both? The 
answer to this question is given by Theorem 3.3: The drift follows a Fick's law in the following sense. The 
configuration space of the n-dimensional ball is Q = SO(n) xR m_1 and there are n—1 constraints corresponding 
to the directions in the table K n_1 . By a symmetry reduction argument (compression) one can eliminate the 
M n_1 -factor and the deterministic motion of the ball can be described by the geodesic equations of the so-called 
non-holonomic connection V nh on SO (n). With respect to this connection one can now show that the process 
on SO(?i) describing the balls stochastic motion is a non-holonomic diffusion separating into a drift- and a 
martingale-term. See Definition 3.2. Theorem 3.3 says that this V nh -drift equals 
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where J\f is the preserved density (3.5) of the deterministic ball and the gradient is computed with respect to 
the kinetic energy metric (3.4) of the ball. 

In particular, when the ball is homogeneous there is no drift. Moreover, it is shown in Corollary 3.4 that the 
homogeneous ball's stochastic process factorizes to a Brownian motion on 'the ultimate reduced configuration 
space' 5 n_1 = (Q /R n_1 ) /SO(n — 1). This is in analogy to the corresponding deterministic case where the 
motion is Hamiltonian at the ultimate reduced level. (See [12].) 

In 3 dimensions the drift — ^grad^ log TV" does not have an angular velocity component about the vertical 
axis of the ball in the space frame. For dimensions n > 3 it turns out that this property is related to the 
Hamiltonization of the deterministic system: If the inertia matrix describing the balls mass distribution satisfies 
the Hamiltonization condition 3.7 then the drift does not have an angular velocity component about the vertical 
axis in the space frame. On the other hand, the drifts angular momentum about the vertical axis is always 0, 
regardless of the dimension or the mass distribution. 

Section 2 starts by collecting some definitions and facts from stochastic differential geometry as presented 
in [14, 10]. Then we rehearse the basics of the stochastic Hamiltonian mechanics and their symmetries as 
introduced in [17, 18]. In Section 2.C these ideas are transferred to describe stochastic G-Chaplygin systems. 
It is noticed that the reduction of symmetries, termed compression in this context, works naturally. 

In Section 3 this construction is applied to the n-dimensional Chaplygin ball. First some facts about the 
deterministic system, such as the preserved measure and Hamiltonization, are recalled. Then the Hamiltonian 
construction of [17] of Brownian motion on the configuration space Q = SO(?i) x R™ _1 is reviewed. Section 3.C 
makes the constraint forces act on the Brownian motion according to the recipe of Section 2.C. Thus a 
constrained stochastic motion is obtained and the above mentioned Theorem 3.3 is found. Finally, we note 
that we can also treat the cases of angular jiggling only or horizontal jiggling only. The latter corresponds to 
the Chaplygin ball sitting on a table which undergoes a translational Brownian motion. The 3-dimensional 
version of this case has been considered in [19]. 

2. Stochastic Chaplygin systems and reduction of symmetries 

2. A. Stochastic geometry. We state some notions from stochastic differential geometry. The references we 
used here are [10, 14]. See also the appendix of [17]. 

Let M be a manifold, let (0, J-, {Ft ■ t > 0}, P) be a filtered probability space, and let Y : R + xfl->M be 
an adapted stochastic process. (We consider only continuous processes.) 

The process T is called a semi-martingale if / o r is a semi-martingale in R for all / G C°°(M). See [10, 
Chapter III]. 

The definition of a martingale in M depends upon a choice of a connection. Let V be a connection in 
TM -> M. Then the Hessian of V is defined by 

Hess v (/)(A, Y) = XY(f) - V x Y(f) 

for X, Y G X(M). This is bilinear in X and Y but not symmetric, in general, since Hess v (/)(X, Y) — 
Hcss v (/)(r,X) = -Tor v (A,r)(/). 

Definition 2.1. A semi- martingale r : R + x O -> M is said to be a V-MARTINGALE if, for any / G C°°(M), 

foTt-foTo- f Ress v {f ){aT s ,aT s )ds 
Jo 

is a local martingale in R. 

In [10, Chapter IV] this is stated in terms of a torsionless connection but it is noted that one can also allow 
for connections with torsion, since it is proved ([10, (3.14)]) that J* Hess v (/)(<ir s , dT s )ds depends only on the 
symmetric part of Hess v . The situation is similar to the notion of a V-geodesic. This also depends only on 
the torsionless part V — ^Tor v of V. 
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Let M and N be manifolds, let (f2, J 7 , {T t ■ t>0}, P) be a filtered probability space, and let X : R + x£l — > N 
be a semi-martingale. A Stratonovich OPERATOR S from TN to TM is a family of linear linear maps 

S (XtV) : T X N — ► TyM 

which depends smoothly on i; 6 ]V and y G M. In other words, <S is a section of T*N <S> TM —> N x M. A 
Stratonovich differential equation for a semi-martingale T : K + x f2 — >• M is written as 

5T = S(T,X)SX. 

See [10, Chapter VII] for the precise meaning of this equation as well as existence and uniqueness (up to 
explosion time) of solutions. 

Assume now that N = R x R n and that X : R + x O -> R x R n , (t,uj) \-> (t, W t (oj)) where W denotes n- 
dimensional Brownian motion. Let Xq, X\, . . . , X n be vectorfields on M and define the Stratonovich operator 

S{x,y) :R n+1 -^T y M, (t, w 1 , . . . , w n ) ^ tX (y) + (tf). 

Let / G C°°(M). Then f oT :R+ xfl^R satisfies 

= df(T).S(X,T)6X = 4f(X).X (T)5t + 'Z4f(X).X i (X)8W i . 

Hence T defines a diffusion in M by [14, Capter V, Thm. 1.2]. The generator of this diffusion is the second 
order differential operator A given by 

Af = J0/ + 5EW 

where / G C°°(M). 

Definition 2.2. If M is equipped with a Riemannian metric /z and A is the associated Laplacian then the 
diffusion is called Brownian motion in (M, /x) if 

A=IA. 

This definition agrees with the one given in [10, (5.16)] or [20, (8.5.18)] but differs slightly from the one in 
[14, Chapter V, Def. 4.2] where it is required that n = dimM. 

2.B. Stochastic Hamiltonian systems. This section presents some of the concepts elaborated in [17, 18]. 

Let again N = R n and consider a Poisson manifold (M, {.,.}) together with an K n -valued Hamiltonian 
function h = (hi, . . . , h n ) : M — > W 1 . Let X : R + x — > R n be a semi-martingale. The associated stochastic 
Hamiltonian system is given by the Stratonovich equation <5r = S(X, T)SX where S is defined in terms of the 
Hamiltonian structure, that is, 

S(x, y) : R n — > T V M, (x\...,x n )^ ^ (vV 

where the Hamiltonian vectorfield X^ is the vectorfield corresponding to the derivation {/, .}. When (M, iS) 
is a symplcctic manifold then one uses the Hamiltonian fields defined by {(Xh^u) = dhi. 

These systems allow for a symmetry reduction analogous to classical mechanics. We state the symplectic 
version of this theorem of [18, Section 6]: 

Theorem 2.3 (Stochastic Hamiltonian reduction). Let (M,uS) be a symplectic manifold with Hamiltonian 
h : M R n and stochastic component X as above. Assume that (M,u>,h) are invariant under the free and 
proper action of a Lie group G such that a coadjoint equivariant momentum map J : M — > g* exists. Fix a 
level A G Q* . 

Then J _1 (A) is invariant under the flow of S. Moreover, S induces a Stratonovich operator S\ from TR n 
to T(J~ 1 (\)/G\) and solutions of S with initial condition in J _1 (A) project to solutions of S\. The induced 
operator is given by 

S\(x,y Q ) : (x 1 ,...,x n ) 1 — > J2 X hf(yo)x % 

where x G R n , yo G J~ 1 (X)/G\ and X h \ is the Hamiltonian vectorfield with respect to the reduced symplectic 
form on J~ 1 (\)/G\ of the induced function . 
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2.C. Stochastic G-Chaplygin systems. A (deterministic) G-Chaplygin system consists of a Riemannian 
configuration manifold (Q, fj,), a Lie group G acting freely and properly by isometries on (Q, fj,), and a horizontal 
space V of the principal bundle 7r : Q -» Q/G. Hence V is the kernel of a connection form A : TQ — > g. The 
Lagrangian of the system is the kinetic energy, i.e., L(q,v) = ^/j, q (v,v). In general, V is not //-orthogonal to 
the vertical space kerT-7r. We will henceforth identify TQ — T*Q via fi. See also Section 4. 

Let Jq ■ TQ — ¥ q* denote the standard equivariant momentum map associated to the lifted G-action on TQ. 
Given a G-invariant function h : TQ — > R one may use Noether's theorem to conclude that the Hamiltonian 
vectorfield Xh is tangent to Jq 1 (0) and, moreover, is pro jectable for JJ^O) -» J G 1 (0)/G = T(Q/G). This 
is why results such as Theorem 2.3 work. Since T> is a perturbed version of Jq (0) we look for a substitute 
construction. 

Consider the horizontal space associated to the pulled back connection t* A = A o Tt : T(TQ) — > TQ g 
of the tangent lifted G-action on TQ, 

T := ker r* A C TTQ. 

By assumption T> is also G-invariant. Thus we can consider the restricted G-action on T> and the associated 
connection l*t* A : TV — > q where l :T> TQ is the inclusion. Define 

C := kertVM C TV 

to be the horizontal space of the principal bundle V -» V/G = T(Q/G) 

According to [1] we can decompose T along V as 

T\V = C © {T\V) n 

where (J-\V) n is the ri-orthogonal of T\V in TTQ\V. In particular, the fiber-wise restriction of to C x C 
is non-degenerate. For z G V define the projection 

(2.1) P Z :T Z {TQ)^F Z ^C Z 

where we first project along the vertical space of the G-action on TQ and then along T z . Moreover, for a 
fc-form <j) on TQ we denote the fiber- wise restriction of L*<p to II fe C by cjf . For a function h : TQ ->Rwe may 
thus define the vectorfield X% on V with values in C by the formula 

(2.2) P z X h {z) = (fl c )-\dhf z =-. X c h (z) 

where z 6 P. The link to to the non-holonomic system introduced above is the following: Via the Legendre 
transform the dynamics of the non-holonomic system (Q,V,L) can be equivalently described by the triple 
(TQ, fi c , H = \n{p,p)) together with equation (2.2). Thus the Lagrange multipliers have been encoded in the 
two-form f2 c or, equivalently, in the projection P : T(TQ)\V — > C. The idea is that a non-holonomic system is 
a Hamiltonian system acted upon by constraint forces. The effect of the forces is described by the projector 
P : T(TQ)\V -> C. 

In addition to this structure consider now a semi-martingale X : R + x £1 — > R™ and a Hamiltonian function 
h = (hi)i : TQ —> R" as above. The associated stochastic non-holonomic system is given by the Stratonovich 
equation ST = S C (X, T)SX where the Stratonovich operator «S C : R™ x V Hom(T]R Ar , C) is given by 

(2.3) S c {n,z) :R" C z , (x\ . . . ,x n ) — >• ^P z X hi [z)x l = 

Thus the non-holonomic Stratonovich operator arises, by applying the constraint forces, as a projection of the 
Hamiltonian Stratonovich operator into C. When the h l are G-invariant we refer to the collection (Q,V,h = 
(hi)i,X) as a stochastic G-Chaplygin system. 

Proposition 2.4 (Compression of stochastic G-Chaplygin systems). Let (Q,V,h = (hi)i,X) be a stochastic 
G-Chaplygin system with £l c as above. Then the Stratonovich operator (2.3) compresses to a Stratonovich 
operator <S nh from TR™ to V/G = T(Q/G) which is given by 

S nh (n, z ) : R" — ► T Za (T(Q/G)), (x 1 , x n ) .— ► J2Kf(^ 

where zq £ T(Q/G) and h® : T(Q/G) — > R is the function induced on the quotient from the invariant function 
L*hi. Moreover, solutions of (2.3) project to solutions of S nh . 
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Proof. Everything is entirely analogous to the proof of [18, Theorem 6.7] with the only difference that now one 
uses Proposition 4.1 instead of the usual symplectic reduction theorem. □ 

We think of ST = S nh (X, T)SX as the equations of motion of the system (Q, T>, h, X). 

Proposition 2.5 (Ito representation). Let f e C°°(T(Q/G)). Then the Ito representation of the equation 
S(foT) =S nh {X,f oT)SX is 

d(f oT) = J2 XfifWdX* + \ X^X^f(T)[dX\dX"} 

i i,j 

where X 1 = pr ! o X . 

Proof. This follows exactly as in the proof [17, Proposition 2.3]. It is only necessary to notice that this proof 
docs not depend on whether or not the non-holonomic bracket Xg h f = {g, /} = — {/, <?} nh satisfies the 
Jacobi identity. The only property of the Poisson bracket which is used in [17, Proposition 2.3] is the Leibniz 
rule and this feature is evidently shared by the non-holonomic bracket. □ 

3. The stochastic Chaplygin ball 

3. A. The deterministic system. For background on the Chaplygin ball we refer to [6, 8, 9, 11, 12, 16]. The 
configuration space of Chaplygin's n-dimensional rolling ball is Q = K x V where K = SO(ro) and V = K n_1 . 
The no-slip constraints are given by the distribution T) = (A + pr 2 ) (0) C TQ where 

A = J2v a ® e a 

where ei,...,e„_i is the standard basis on V and we stick to the following conventions: 1 TK = K X t 
is trivialized via left-multiplication and t is equipped with the Ad-invariant inner product (.,.}. 2 Let H = 
SO(n — 1) C K be the stabilizer in K of the n-th standard vector e„ such that H acts on V in the natural 
way. We decompose t = f) © f) with respect to the Ad-invariant inner product (., .) on {. With respect to this 
inner product wc introduce an orthonormal system 

Yq, a = 1 , . . . , k = dim f) and Z a , a = 1 , . . . , n — 1 

on t such that Y a £ f) and Z a G f) . Associated to this basis we define the right invariant vector fields £ a and 
C a . In the left trivialization these read 

= Ad( S " 1 )y Q and Us) = kd{s- x )Z a . 
Dually we introduce the corresponding right invariant coframe 

P Q H£ Q ,.} and »j = (&,.>. 

The Lagrangian is the function 

L :TQ = K xtxTV — > R, (s,u,x,x') i — > \{lu,u) + \{x\x'). 

where I is the inertia matrix in body coordinates. The rolling ball with the no-slip constraint is the non- 
holonomic system described by the data (Q,V,L) where the equations of motion follow from the Lagrange- 
dAlembert principle. However, we will not have much use for the Lagrange function below since we will only 
perturb the resting ball. Note also that we overload the symbol (., .) by using it for the Euclidean inner product 
on V as well as for the Ad-invariant structure on t. 

From a structural point of view the decisive feature of the Chaplygin ball is that its constraints are given 
by a connection pr 2 + A : TQ — > V on the (trivial) principal bundle V K X V -» K where V acts on 
itself by addition. Thus T> is the horizontal space of this connection. However, T> is not /i-orthogonal to 
the vertical space of the bundle. The fact that the system is non-holonomic is reflected in the non-flatness, 
Curv" 4 = dA ^ 0. 



In Section 2.C it was the connection form A + pr 2 which was called A.. 

The space Xl (K) of left invariant vectorfields and the Lie algebra t will also be identified without further notice. 
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Compression of the Chaplygin ball system yields the almost Hamiltonian system (TK, f2 n h, TL C ) where f2 n h 
is described in Proposition 4.1, TK and T*K are identified via the induced metric 

(3.4) fi Q ( Ul ,u 2 ) s = (lu 1 ,u 2 ) + (A s (ui),A s (u 2 )) = ((I + A* s A s )ui,u 2 ), 

and TL C = i//o(it, u) s . The metric is the sum of a left invariant and a right invariant term. Thus it constitutes 
an L + i?-system, see [15]. Note the useful formula A*A{u) = XXCa,w)Ca- 

The compressed system (TK , Sl n h, H c ) is further invariant under the lift of the left multiplication action 
of H on K. Physically this corresponds to rotation of the ball about the e„-axis in the space frame. This 
is an inner symmetry and gives, by the non-holonomic Noether theorem, rise to a conserved quantity. This 
quantity is just the standard momentum map Jh ■ TK — > t)* = t), (s, u) i-4 ^2p^(lu)Y a of the (co-)tangent 
lifted iJ-action on TK. 

The Chaplygin ball shares an important feature with Hamiltonian systems. Namely, it possesses a preserved 
measure ([6, 11]). At the compressed level -the TK -level- the density Af : K — >• R of this measure with respect 
to the Liouville volume on TK is 

(3.5) Af(s) = (det/x (s))~*- 

This function plays the central role in all questions of Hamiltonization of the system. Note that J\f is H- 
invariant and thus descends to a function J\f : K/H = 5 n_1 — > K. For further reference we also record 
that 

(3.6) rf(log^) = E([i«o" 1 Ca J Ca] J a)ry 6 . 

In [12] it is proved that f2 n h can be replaced by 17 = fl K — gXXtCai Cb], -)v a A if without altering the 
equations of motion, that is, i(X^ h )f2 n h = i(X^)£l = dH c . Now the new system (TK,n,H c ) has the same 
dynamics but has the advantage of being liable to reduction with respect to the internal symmetry group H : 
Let Jh ■ TK — > f)* be the momentum map introduced above, A 6 fj* and t : J^ X (A) <—t TK the inclusion. 
Then l*Q descends to an almost symplectic two form fl\ on J [{ 1 (X)/H\. (t*O n h is not horizontal for the 
projection onto JJ 1 1 (X)/H\.) In this way one can recover the Hamiltonization of the 3-dimensional ball of 
[4, 5]: It is shown in [12, Proposition 4.4] that d(Aftl\) = if n = 3. Thus, for n = 3, the rescaled vectorfield 
Af~ x X^ is Hamiltonian with respect to (J H 1 {X)/H X = TS 2 , A/T2 A , U c ). Moreover, the homo gencous ball, 
I = 1, is Hamiltonian at the J^ 1 (X) / H\-\evel for any dimension n. It is interesting to notice that none of these 
statements hold at the Tif-level. 

In [13, Section 4. A] the following Hamiltonization condition is proved for arbitrary dimension n: Let A = 
£ f)*. Then d(AfCt Q ) = if and only if 

(3-7) (n - 2)(^ 1 Q, [C b ,Cc]) = Ea^o'Ca, KbXa]S c , d - [Cc,Ca]M 

which is an algebraic condition on I. For n = 3 this condition is trivially satisfied. In the stochastic context 
this condition appears in Theorem 3.3 below. 

3.B. Brownian motion on the configuration space. We follow [17] to construct Brownian motion on 
Q = K x V. Let V M be the Levi-Civita connection of 

fi({u,x'),{v,y')) = (Iu,v) + (x',y'). 

Thus for iijWfif we have 

V£« - (Vi«,0) = [\[u,v] + ^([uM + [«,!«]), 0) 

where V 1 is the Levi-Civita connection of the left-invariant metric defined by (I., .) on K. Note that we identify 
6 = t* via (., .} whence I : I -> t* = I. Let 

Vi, . . . ,v m , m = dimt = ^n(n — 1) 

denote a basis which is orthonormal for (I., .). Thus Vi, . . . , v m , ei, . . . , e n _i is a left invariant frame on Q 
which is orthonormal with respect to the left invariant metric /i. (Remember that we identify TQ and T*Q 
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via /x.) Let the functions Hq, Hi,F a : TQ — > R be given by 

(3.8) H (s,u,x,x') = —3 (In, Hi(s,u,x,x') = (Iu,Vi), and F (s,«,x,x') = {x ,e a ). 
Consider the semi-martingale 

I:K + x!]^lx R m x R"" 1 , (i.w) .— > (t, B}(w), . . . , -Bf(w), W^w), . . . , W t " _1 M) 

where B l ,W a are m + n— 1 independent Brownian motions. The Stratonovich stochastic differential equation 
which is associated to these data is 

(3.9) sr = s Ham {x,r)sx 

where the Stratonovich operator from T(R x R m x R n_1 ) to TK x TV is dchned by 

5 Ham (n; s, u> x, x')(t, b 1 , . . . , b m , w 1 , . . . , w"" 1 ) = X Ho (s, u, x, x')t + J2 x h, (*, u, x, x')b' + J2 x F a (*, u, s, 

with X// denoting the canonical Hamiltonian vectorfield of a function H : TQ — > R. Using the Ito repre- 
sentation of this equation [17] show that the solutions F project via r : TQ — > Q onto Brownian motion on 
Q. 

Using again the setting of [17] and Theorem 2.3 it is easy to see the following. Consider the U-action on Q 
as above. Let Jy 1 (0) = {(s, u,x,0)} be the 0-level set of the standard momentum map Jy ■ TQ — > V* = V 
of the lifted V-action on TQ. Then the Stratonovich equation (3.9) induces a Stratonovich equation (JITo = 
So(X, Tq)SX on Jy 1 (Q)/V = TK and solutions T with initial condition in J^ 1 (0) project onto solutions To. 
Moreover, tk o To is a Brownian motion on K where tk ■ TK — > K. 

When we regard T> as perturbed version of Jy 1 (0) we can ask how much of this observation remains true? 
This is the content of Section 3.C. 

3.C. Constrained Brownian motion and compression. Wc now force the Brownian motion on Q to 
satisfy the constraints induced by T>. In accordance with Section 2.C we do so by applying the constraint 
forces to the Stratonovich operator from (3.9). Thus we are concerned with the equation 5T = S nh (X,T)SX 
where X and Ho, Hi, F a are as above and 

S c (n, z) = P z S Ham {n, z) : R x R m x R"" 1 — ■+ C z C T Z V 

where P was defined in (2.1) and z = (s,u,x, —A s (u)) G T>. The functions Ho, Hi, F a are U-invariant and 
compress to functions ho, hi, f a given by 

h (s,u) = H Q (s,u,x, -A s {u)) = -|/i((-u, 

= -|(i«,Ev*3) - -ImoK/^iEv^), 

hi(s,u) = (lu,Vi) =/io(u,/io 1 Ivt) ) 

fa(s,u) = -(A s (u),e a ) = -H (u,LJ,o 1 Ca(s)). 

Notice that ho and hi are left invariant while the f a are right invariant. The compressed non-holonomic 
Stratonovich operator is now of the form 

<S nh (n; s, u)(t, b 1 ,..., b m , w 1 ,..., w^ 1 ) = X^(s, u)t + £*/" h (s, u)b l + £X£ h ( S , u)w a . 

We think of solutions of ST = S nh (X, T)8X as non-holonomic diffusions. This is in analogy to [2, Chapter V] 
where Hamiltonian diffusions are considered in a similar manner. 

For a function / £ C°° (K) , viewed as a function on TK via pull-back, we have 

(3.10) X^(f) = -\df.iXo X lTyl^ X^{f) = df.^Ivi, and Xj*(f) = -df.^Ha- 

Let r be the solution semi-martingale to ST = S nh (X, T)SX and let r : TK — > K be the projection. Then 
r o r solves 

S(t o F) = Tt.ST = -^o ll E^lM T ° r )^ + X>o 1] M T ° r )^ 1 - EMo X Ca(r ° T)5W a . 
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According to [14, Chapter V, Theorem 1.2] this means that the semi-martingale r o T defines a diffusion in K 
whose generator is the second order differential operator 

(3.11) - l^nEvJ^ + ±E(mo ^Xmo ^) + iEO^OO^Ca). 

To identify the drift of the diffusion roT a connection is needed. We introduce the non-holonomic connection 
which is explained in [7, Section 5.1.1] : 3 Let Pr M : TQ — > V denote the projection onto T> along the /x- 
orthogonal V 1 of V. Note that ^ kcrT?r where n : Q -» K is the projection. Let hi" 4 : X(K) -4 £(Q; 2?), 
hi (s, a:) (it) = (s, u, cc, — _4 s (w)) be the horizontal lift map associated to A. Given A, Y G X(A) the non- 
holonomic connection is prescribed by 

(3.12) VfY = TTrPi tl V^ Ax {h\ A Y). 

This connection is metric, i.e., V nh fio = 0, and its geodesic equations are exactly the equations of motion of 
the non-holonomic system described by (TK, f2 n hj He)- However, in general, V nh will have non-trivial torsion. 

Lemma 3.1. For ii,wet = Xl(K) we have 

where V^f = §[it, v] + ^I^ 1 ([u,Iv] + [w,Iu]). Its torsion is given by Tor nh (w, v) = (1q 1 A*A[u,v\. 

Note that V" h w is not left invariant any more. At the compressed level the equations of motion of the 
Chaplygin ball write as u' + /i " 1 [u,Iu] = 0. In 3D this corresponds to [8, Equation (3.5)]. 

Proof. For u,v <E t we need to compute 

P ^K^-Au)( v > ~M = Pr M (V^ - u.Av) = Pr"(Vj,« - A[u, v}). 
Now note that (w, X) G T>^ if and only if w = I^ 1 A*X. We have to solve 

= A(W J u v + I- 1 A*X) - A[u, v] + X 
for I~ 1 A*X. The solution is found to be given by fi I- 1 A*X = A*A{[u,v] - V l u v). Therefore, 

(3.13) Vf v = W J u v + ^ 1 A*A([u, v] - W l u v) = ^ 1 A*A[u, v] + (i^Wlv. 

□ 

Recall that the Hessian of V nh is defined by 

Hess nh (/)(A,r) = XY(f) - VfY(f) 

for X, Y 6 X{K) and / G C°°{K). 

Let / G C°°(K) and ST = S nh (X,T)SX. By Proposition 2.5, equations (3.10) and [dB\dB*] = S^dt, 
[dW a ,dW b ] = 6 a > b dt, as well as [dB\ dW a ] = we have the Ito equation 

d(f or) = ( - ^lYylxn + |e(mo x H)(mo + lE^O^C.))/^* 
+ E(^o lfo 0/(r)rfi? 1 - EO^CJ/OW 
= ( - W l Yyl^ + IEV^i^^ 1 !^) + !EvJ lf >o x Ca) 

+ iEHess nh (Mo 1 fei ! /»o 1 ^) + lEHess^o'Ca.Mo'Ca))/^ 
+ ElftV)/^ - EO^O/OW, 

which also confirms (3.11). Having split the generator into first and purely second order part it makes sense 
to say what we mean by drift. 

Definition 3.2. The vectorficld 

-K'nEv^ + lEvJ^O^ifi) + |EV^ 1Co (mo *&) 

is called the drift of the diffusion toT with respect to V nh . 

■^Contrary to [7] we only use the projected version of the non-holonomic connection. 
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Note that t o T is a V nh -martingale if and only if the V nh -drift vanishes. See also [10, Theorem (7.31)]. 

Theorem 3.3. Let T be a solution of the Stratonovich equation 8T = S nh (X, T)8X and let t : Tlx — > K be 
the projection. 

(1) Then, with respect to the non-holonomic connection V nh introduced in (3.12) ; the semi-martingale r o T 
defines a diffusion on K whose drift is the gradient — igrad W) (log J\f) where Af is the density function 
defined in (3.5). 

(2) The drift — ^grad' 10 (log J\f) is horizontal with respect to the mechanical connection, Hor moc = (keiTu)^ - 1 , 
on the principal bundle k : K -» K/H = S n . Ifl satisfies the Hamiltonization condition (3.7) then the 
drift is also horizontal with respect to the principal bundle connection E/ 9 " C2> Y a : TK — >• f) on k : K -» 
K/H. 

Item (2) means that the drift's component of angular momentum about the vertical axis in the space frame 
vanishes, and when the Hamiltonization condition holds then the same is true for the component of angular 
velocity about the vertical axis. For n = 3 this condition is always satisfied. 

Corollary 3.4. Let T be a solution of the Stratonovich equation 5T = S ah (X, T)SX and let r : TK —> K and 

k : K —» K/H = S'""^ 1 be the obvious projections. Suppose that 1=1 (i.e., the ball is homogeneous). 

(1) Then roT defines a martingale in K with respect to the non-holonomic connection. 

(2) The process kotoT is a Brownian motion on if?™ -1 whose generator is \ times the Laplacian of v , where 
v is the metric on induced from the left H -invariant metric hq = ((1 + „4*.4)., .) on K . 

Note that the restriction of /io to Hor = Hor mcch = span{£ a } equals twice the restriction of the biinvariant 
metric. This corollary is intuitive but nevertheless not obvious since the dynamics at the compressed level can 
never be described by a Hamiltonian reduction procedure. This is because A is not the mechanical connection, 
even if the ball is homogeneous. (Compare with [12, Corollary 4.3].) Thus it does not fall in the category of 
[17, 18]. 

For reference we note the formula 

(3.14) V.fi (u) = V.J2(Ca,u)Ca = £(Ca, [v, U ))( a - E(Ca, u)[v, Ca] = A*A[v, u] - [v,A*Au]. 

Proof of Theorem 3.3. Let / £ C°°(K) which we regard via pull-back as a function on TK. According to 
Definition 3.2 we need to show that 

(3.15) - Mo^lEV 1 ^ + E^ lVi (^i) + EVjx c >o \a) = -grad^QogAO- 
Claim: 

(3.16) EVj^Jft 1 ^) = Mo 'lEV^i - grad^(logAT). 
Indeed, we use (3.14) and the fact that u = J2(^ v i> u ) v i to see that 

EVj lhi ( ft - l bi) = £<i/io 1 fo,,^)(i A1 o Hvuv^v^vk + E^oV^iKMWV.^H 

= E Mo 1 [/^O 1 Iv i > ^0 1 lv i ] + E 1 lv i ' v j ) (^0 1 [ V 3 J A* A(1q 1 1 V t ] , V k )v k 

= EMo Vo ^jJMo ±Iv i} + EMo Vo 1 M,A* A^lVi] 

= EMo Vo llv i> Iv i) 

= EMo 1 [ v iJvi] - J2t l a 1 [l l o 1 A*Avi,lvi\ 

where we use /Iq 1 I = 1 — /-i~ 1 A*A in the last equation. Notice that hv,,!?^] = IVj,.?^. For the gradient part of 
claim (3.16) we consider 

= E([Mo 1 Ca,Ca],Cfc) =<l0g^)C6 
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by (3.6). Similarly it is true that (EI^o 1 A*Avi, Ca) = E([Mo CaL C<*)- Using the property [Z a ,Y a ] = 
S a ,bZ c — § a ,cZb with Y" Q = [Zb, Z c ] it is easy to see that 

Thus ^^o 1 !^ 1 - 4 *- 4 ^. 1 ";! = /"o 1 E(#°g^)C6)C6 = grad Mo (logAf) and claim (3.16) follows. 
Claim: 

(3.17) V^ 1Ca ( Mo - 1 Ca) = 0. 

Again we use (3.14) to see that 

V^Gtio'Ca) = Mo'lMo'CaJ^o'Ca] +E(H>^ 1 Ca)(^.(fcfc^0 1 Ca))^ 

= MO 1 \M0 1 Ca , IMO 1 Ca] + E : M(7 1 Ca ) (If* , Mo 1 . -4* -^^0 * Ca] K 
- J] (!«,• , ^ 1 Ca ) (lUjfe , Mo 1 > Ca] }Vk 

= Mo 1 [/ U 1 Ca,lMo 1 Cj +/i(7 1 [M(7 1 Ca,^*^Mo 1 Ca] - Mo 1 [Mo ^a, Ca] 
= 0. 

Now (3.16) and (3.17) imply (3.15) which shows part (1) of the assertion. 
For (2) one checks that (3.7) yields, for £ a = [£ c , Q], 

(grad" (log AO, Ca) = EW'lCcCd], ([Mo 1 Ca,Ca],C6)Cb) 

= E(M(7 1 [Cc,Cd],C b )(M : 1 Ca,[Ca,C o ]) 

= (MoT'Cd, [Co, Cd])<Mo 'Co, [Cc Ca])(n - 2) + (/i^Cc, [Co, Cd] ) (m ~ 1 Cd , [Cd, Cc]>(n - 2) 
= 

where we also have used that E([Mo Ca, Ca], Ca) =0. □ 

Proof of Corollary 3.4- Part (1) is clear. 

Concerning part (2) let / G C°°(S'™ _1 ). According to Definition 2.2 we should to show that the generator 
(3.11) satisfies 

(iEGCOGuoV) + |E(mo 'CcOOuo x Ca)K/ = |«*A"/ 

where A 1 ' is the Laplacian associated to Indeed, we find tjE(Mo Viii^o v i) K * f + !>E(Mo 1 Ca)(/ i o" 1 Ca) K */ = 
flECaCa^*/. Now {-^=Ca} is a horizontal orthonormal frame for /zo | (Hor x Hor) where Hor is the /x -orthogonal 
to kerTre. Therefore, 

lECaCaK*/ = iTr hOT Hess«'( K *J) = |«*A"/ 

where Tr hor denotes the trace computed with respect to horizontal fields only and Hess Mo is the Hessian of the 
Levi-Civita connection on (K, /j, ). The equation ( a ( a K* f = Hess Mo («*/)(Ca, Ca) is justified by the observation 
that V^°Ca = where V* is the Levi-Civita connection of the right invariant metric /io- (In fact, the restriction 
of V nh to Hor x Hor equals the restriction of V W) . ) □ 

In the homogeneous case the above construction yields a Brownian motion on S" -1 in a manner similar to 
the one described in [14, Chapter V] by the notion of rolling the sphere S*™ -1 along a Brownian motion in 
by means of the Levi-Civita connection. The difference is that [14] start from Brownian motion in W n ~ 1 while 
we started from Brownian motion in R m x with m = dim so (n) — "^" 2 ~ 1 - > . One can recover the setting of 

[14] by setting (Hq, Hi, ... , H m ) — in (3.8). Then, with I = 1, we obtain a diffusion kot oT which is driven 
by Brownian motion (W 1 , . . . , W"^ 1 ) in R™ -1 and the generator of which is given by 

^E(Mo 1 Ca)(Mo" 1 Ca) K */=^*AV 

where / £ C°° (S""" 1 ). Referring to the interpretation stated in the introduction this means that the Chaplygin 
ball is subjected to horizontal jiggling but there is no angular jiggling. Equivalently, the ball sits on a table 
which undergoes a translational Brownian motion. Compare [19]. Alternatively, we can set (fi, . . . , F n — i) = 
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in (3.8). Then there is only angular jiggling and the diffusion roT is driven by (B 1 , . . . , B m ). By (3.17) the 
drift remains the same as in Theorem 3.3. 

It seems that the notion of a stochastic non-holonomic system has been hardly investigated in the literature. 
We finish by asking the following questions. 

(1) Does an analog of Theorem 3.3 hold for general G-Chaplygin systems when there is a preserved measure? 
What can be said about the drift if there is no preserved measure? 

(2) Which is the precise relationship between G-Chaplygin systems with preserved measures and measure 
preserving 'Chaplygin diffusions'? Preservation of measure by diffusions is studied in [14, Chapter V]. 

(3) Is there a time change or Girsanov type argument to eliminate the drift in Theorem 3.3 or to make it even 
a Brownian motion? Is this related to the Hamiltonization of the deterministic problem? 

4. Appendix: G-Chaplygin systems and symmetry reduction 

The purpose of this appendix is to shortly introduce and motivate the notion of a G-Chaplygin system and 
to state Proposition 4.1 which explains the symmetry reduction of such systems. This reduction is termed 
compression ([9]) to distinguish it from symplectic reduction. These concepts are closely related and compres- 
sion can be viewed as a perturbed version of its symplectic counterpart. At the same time, however, there 
are fundamental differences; symmetries behave differently in non-holonomic mechanics and do not necessarily 
give rise to conserved quantities, and there need not exist a preserved measure ([7, Section 5.4]); all this is 
related to the question of closedness of the form f2 n h defined in Proposition 4.1. See [1, 3, 7, 9, 12]. 

A non-holonomic system is a triple (Q, 2?, L) where Q is a configuration manifold, L : TQ I is a 
Lagrangian, and T> C TQ is a smooth non-intcgrable distribution which is supposed to be of constant rank. 
The equations of motion for a curve q(t) which should satisfy q' £ T> are then stated in terms of the Lagrange 
d'Alembert principle. Suppose there is a Riemannian metric fi on Q such that we have an isomorphism 
TQ = T*Q and assume that L is the kinetic energy Lagrangian. In this case there is also an (almost) 
Hamiltonian version: continue to use the symbol /i to denote the co-metric and consider the Hamiltonian 
T-L(q,p) given by the Legendre transform of L. Since T> is of constant rank there is a family of independent 
one-forms <j> a 6 £l(Q) such that T> is the joint kernel of these. In terms of coordinates (q l ,Pi) the equations of 
motion are 

where the A Q are the Lagrange multipliers to be determined from the supplementary condition that p(4> a ,p) = 0. 
With := (q',p') we may thus rephrase the equations as 

(4.18) i(X%t)n = <M + Y.^aT*4> a 

where fl = —d6 is the canonical symplectic form on T*Q and t : T*Q —> Q is the footpoint projection. (The 
notation will become clear below.) 

Let G be a Lie group that acts freely, properly and by isometries on the Riemannian manifold (Q, (i). A G- 
Chaplygin system is a non-holonomic system (Q,L = 1 1 1 • 1 V) that has the property that I? is a principal 
connection on the principal bundle Q -» Q/G. Thus T> is the kernel of a connection form A : TQ — > g. Notice 
that we do not require A to be the mechanical connection associated to [i. 

Consider C and il c as defined in Section 2.C. Since is, by construction, tangent to M. and takes values 
in C one may now rewrite the equations of motion (4.18) in the appealing format 

i(X%t )il c = (dtff 

where {dT-L) is the restriction of L*dH. to C with i :T> TQ being the inclusion. 

Let /Lto denote the induced metric on S := Q/G that makes it : Q —» S a Riemannian submersion. Identify 
tangent and cotangent space of Q and 5* via their respective metrics. Consider the orbit projection map 



p := Tn\V : V -» V/G = TS. 
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We may also associate a fiber-wise inverse to this mapping which is given by the horizontal lift mapping hi 
associated to A. As already noted in Section 2.C, l*t* A : TV — > q defines a principal bundle connection for p, 
whose horizontal space is given by C. 

Proposition 4.1 (Compression). The following are true. 

(1) il c descends to a non- degenerate two-form r? n h on TS . 

(2) n h = £ls — (Jg ° fil i TgCurv" 4 }. Here fls = —dQs is the canonical form on TS, Jg is the momentum 
map of the tangent lifted G-action on TQ, Curv" 4 G fl 2 (S, q) is the curvature form of A, and t$ ■ TS — > S 
is the projection. 

(3) Let h : TQ — > K be G-invariant. Then the vectorfield 

xi-.= (tfrw, 

where (dh) is the restriction of t*dh to C, is p-related to the vector field X£ h on T*S defined by 

i(x£)n nh = dh 

where the compressed Hamiltonian, ho ■ T* S = TS — > R is defined by ho := h o i o hi , with hi" 4 denoting 
the horizontal lift mapping. 

In general, J7 n h is an almost symplcctic form, that is, it is non-degenerate and non-closed. See [1, 9, 12]. 
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